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. Abstract. A generalization to the Kraichnan gaussian, white noise in time model of turbulent 

velocity field is proposed. The generalization is designed to take into account the effects of 
stratification in the atmospheric boundary layer by introducing different scaling exponents for 
horizontal and vertical velocity components, leading to strong vertical anisotropy. The inevitable 
compressibility of the model leads one to take into account also the density statistics of the fluid, 
which will immediately lead to a boundary layer like structure in the flow statistics and non zero 
time correlations. Possible applications for passive scalar, tracer and inertial particle statistics 
C"| ' and their anomalous scaling behavior will be discussed. 
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1. Introduction 

Compared to homogeneous, isotropic turbulence, perhaps the most significant difference in 
a stably stratified atmospheric boundary layer is the strong vertical anisotropy manifesting 
as a oc /c -3 mesoscale spectrum (corresponding to oc r 2 spatial behavior in the structure 
function) of the vertical velocity components, as opposed to the usual Kolmogorov-Obukhov 
oc fc -5 / 3 spectrum of the horizontal components (see e.g. [;2J and references therein). The 
vertical velocities are "smoothened" by the effect of gravity pulling the atmosphere toward the 
ground. Recently, much understanding on the nature of intermittency and anomalous scaling 
(or multiscaling) has been acquired by the so called Kraichnan model, where one assumes the 
turbulent velocity statistics to be known and completely determined by a gaussian, white in 
time random velocity field (see e.g. [3] for a review). The purpose of the present work is to 
generalize the Kraichnan model to mimic the behavior of stratified flows. The ultimate goal of 
the project would be to determine how the stratification affects the anomalous scaling exponents 
of the passive scalar, and to study possible clustering of inertial particles. 

2. Kraichnan model 

The Kraichnan model is a gaussian, mean zero, white in time random field which can be 
completely defined via the pair correlation function 

x)i/„(t',x')> = *(At)iV(Ax), (1) 

where S (At) = S(t — t') is the Dirac delta function and D^ u is an incompressible tensor field, 
which can be defined e.g. as 

2V(r) - D Jd P (p2 + L _ 2)3 /2 +5 /2 e • ( 2 ) 



The constant L is the integral scale of the flow and < £ < 2 describes the spatial roughness 
of the velocity field and appears as a scling exponent in the velocity structure function in the 
limit of infinite L as 

<|v(t,x) - v(t',x')| 2 ) oc <5(Af)|Ax| ? . (3) 

One should note that although the infinite L limit of the structure function is finite, the tensor 
field of d2|) diverges as L^. 



The passive scalar evolves according to the equation 



d t 9 - kA6 +v-V9 = f, (4) 

where / = / (t, r) is usually defined as a gaussian, white in time forcing, defined as 

(f(t,r)f(t',r')) = 5(At)C(Ar/L f ) (5) 

with a forcing scale Lf. In the case of an isotropic, large scale forcing, i.e. C — >• C(0) as Lf — > oo, 
the passive scalar equal time structure functions scale as 

([6(t, x + r) - d(t, x)] W ) oc r f (2-0-Aiv ( 6 ) 

with Ajv = N <2(d+2) £ ^(£ 2 ) f° r sma U £• The simple gaussian model of turbulence has therefore 
lead to an intermittency effect manifesting as multiscaling of the passive scalar structure 
functions. 



3. Modified Kraichnan model with stratification 

The idea of the present work is to study the effect of stratification in the context of a Kraichnan 
like turbulence model. For example, one would like to determine if and how the stratification 
affects the anomalous scaling exponents Ajv above. We define a stratified generalization of the 
Kraichnan model as 

for the non- vertical components with i £ 1,2, \i € 1,2,3 and r% is the vertical component. The 
vertical component is taken as 

n^n^J^ —^f^ ^. ( s) 

Here < £ < i] < 2 are the horizontal and vertical scaling exponents, a < 1 is the aspect ratio 
of vertical to horizontal integral scales and the term L^~ v term in the prefactor ensures correct 
dimensionality. In an atmospheric boundary layer one would have rj ~ 2, but we will keep r/ as 
a free parameter. Also included is an ultraviolet cutoff term oc e~ lp , where the length scale I 
corresponds to a viscous dissipative scale. It is often omitted in the Kraichnan model, since one 
is usually interested in scales larger than the viscous dissipation and since the model is finite as 
/ is taken to zero. Here the viscous scale / must be kept finite (for now) because the model is 
not incompressible due to the different horizontal and vertical scaling. Specifically, the quantity 
((V • v) 2 ) diver ges as I approaches zero. The diffusivity ratio behaves as D^/Da oc a 71 , which 
implies vertically inhibited diffusion for small aspect ratios, as is known to occur in real stratified 



turbulence (see e.g. [1]). 



Finite compressibility in the model forces one to reconsider the role of density in the model, 
i.e. it is no longer possible to simply assume a constant density. Instead we take the density as 
another random field beside the velocity field. This is most conveniently achieved by rewriting 
the mass continuity equation 

dtp + V • (vp) = (9) 

in terms of the logarithm of the density (henceforth called the log-density) x = l°g (p/po)> where 
po is some reference density. The reason can be seen from the passive scalar equation, which in 
the case of nonconstant density reads 

d t (p6) + V • (yp9) - kV • (pV6) = 0. (10) 

As usual, the above equations should be interpreted as a stochastic PDE of Stratonovich (mid 
point) type [5]. Taking this carefully into account and going from Stratonovich to ltd prescription 
(see e.g. the appendix A of [I]), using the continuity equation Q and finally taking the limit of 
vanishing molecular diffusivity k, we obtain 

a t e + + £V(o)M d^e - ±£> e = o, (11) 

with the notation Dq = D^^d^dy. The passive scalar is therefore advected by an effective 
velocity field composed of the Kraichnan velocity field v and a gradient of the log-density x- 
The equation for the log-density follows by similar use of the rules of stochastic calculus from 
eq. ©, giving 

dtX + V • v + v • Vx - \d oX + 7 = 0, (12) 
with notation 7 = ((V • v) 2 ) = —D^ jflL ,(0) (symbols after a comma denote derivatives). 



Using the above SPDE, we write the steady state equations for the mean log-density 
g(z) = (x) (which we assume to depend only on the vertical coordinate z = X3), 

Dog = 27, (13) 
which can be readily solved with suitable boundary conditions, yielding 

(x> = g(z) = ( uz ~ \ ) ' ( 14 ) 

where uj 2 = 2((V • v) 2 )/(v 2 ) oc a V (l/L)^l 2 . The solution will therefore tend to zero as L 
approaches infinity. It is perhaps more instructive to consider the mean density rather than the 
mean log-density. The solution for mean density is 

(p) = po cos (uz) , (15) 

which leads one to consider the "height" of the boundary layer as h = ir/2oj oc L^/ 2 . The gradi- 
ent of the mean log-density (|14[) implies that the mean effective velocity has a nonzero vertical 
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Figure 1. Plots of mean density (solid), mean log-density (small dashes) and vertical mean 
effective velocity (large dashes). 



component, pointing towards the mean height. This implies that passive quantities will exhibit 
mean values concentrated at height h/2. The situation is illustrated in Fig. ([1]). Clearly this 
is not a completely realistic density profile, as one might expect since the velocity correlation 
function does not exhibit any height dependence in this approximation. 



The equation for the flux like correlator J-^ (t — t', x — x') = (x(t, x)v^(t', x')) is 

dtT» = S(t - t') [D mv + g'{z)D^] + ~5 T„. (16) 

Since it is the gradient of the log-density that appears e.g. in the passive scalar equation, we 
take a derivative of the solution, which at leading order as L — > oo is approximately 

(d^ X (t, *)vi(t', x')) = -9(At)e^ Atdo d^ [D 3l g'(z) + D iv>v ] (17) 
« 6(At) [c^ + cl v (r/Lf + c%At/L 2 ^] , (18) 

where the c l are some anisotropic unit tensors and r = |Ax|. If we further assume that 
At/L 2 ~£ S> 1, we obtain (x v 'd ^ &(At)(At)~( 5 ~^/ 2 , i.e. the long time time correlation is a 
power law. Similar result for the vertical velocity component can be derived but with £ replaced 
by 7]. 



For the equal time two point function Go(x, x') = (x(t, x)x(£, x')) we have 

( d A + 2d A + Go - d^dlGo - D„ - 7 (g(z) + g{z')) = 0. (19) 

At leading order the last term vanishes. We can then assume homogeneity, which results also in 
the vanishing of the first term, yielding the leading order solution (again with gradients of x) 



(0„x(t,x)$,x(t,xO) w-fyM^A^Cx-x') (20) 
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cLr~ 2 , (21) 



with another anisotropic unit tensor and notation d = d^d^dy = (1)^(0) — D^ u ) d^dy. We've 
also assumed r 2> /, hence the divergence at r = 0. 

The equation for steady state unequal time two point function G(t — i',x, x') = 
(x(*,x)x(t',x')) is in turn, 

d t G = ^D G- 7 g(z'), t>t>. (22) 

The solution, again at leading order, is 

(^x^xRx^x)} ~ -e^^d^Go (23) 

« C % L ^JAJ\ ( 24 ) 

with another anisotropic unit tensor. 
4. Concluding remarks 

The next step in the process would naturally be to consider e.g. the passive scalar and tracer or 
inertial particles in a flow defined by the above model. From the approximative expressions 
above, one can already gain some insight into some properties of these passive quantities. 
First, since there is now a mean effective velocity field pointing toward the height h/2, this 
is where one may expect the particles and passive scalar fields to be concentrated. Second, it 
was demonstrated that the log-density gradient field exhibits a time-like correlation, as opposed 
to the white noise in time correlation of the Kraichnan model. As is well known, the phenomenon 
of clustering of inertial particles requires a non zero correlation time on order for the particles 
to have time to traverse to the edges of eddies. It is therefore reasonable to expect that such 
clustering may occur for this model, simply due to the anisotropic density fluctuations. Also, a 
common critique toward extensions of the Kraichnan model to models with time like correlations 
is that they violate Galilean invariance. Here such a problem does not arise since the log-density 
field is simply determined by the mass continuity equation, which is of course Galilean invariant 
by construction. 
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